EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-EP/98-047 
23 March 1998 



SEARCH FOR A NEW GAUGE BOSON IN tt° DECAYS 

The NOMAD Collaboration 



Abstract 

A search was made for a new light gauge boson X which might be produced in 
7T° — > 7 + A decay from neutral pions generated by 450-GeV protons in the CERN SPS 
neutrino target. The X's would penetrate the downstream shielding and be observed in 
the NOMAD detector via the Primakoff effect, in the process of X — > ir° conversion in 
the external Coulomb field of a nucleus. With 1.45 x 10 18 protons on target, 20 candidate 
events with energy between 8 and 140 GeV were found from the analysis of neutrino 
data. This number is in agreement with the expectation of 18.1±2.8 background events 
from standard neutrino processes. A new 90 % C.L. upper limit on the branching ratio 
Br(TT° —> 7 + X) < (3.3 to 1.9)xl0~ 5 for X masses ranging from to 120 MeV/c 2 is 
obtained. 
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1 Introduction 



Many extensions of the Standard Model such as GUTs fl|], super-symmetric [||, super-string mod- 
els |H and models including a new long-range interaction, i.e. the fifth force (4j, predict an extra 
U (1) factor and therefore the existence of a new gauge boson X corresponding to this new group. 
The predictions for the mass of the X boson are not very firm and it could be light enough 
(Mx <C Mz) for searches at low energies. 

If the mass Mx is of the order of the pion mass, an effective search could be conducted for 
this new vector boson in the radiative decays of neutral pseudoscalar mesons P — > 7 + X, where 
P = 7r°,r/, or 7/, because the decay rate of P — > 7 + any new particles with spin or j proves 
to be negligibly small Therefore, a positive result in the direct search for these decay modes 
could be interpreted unambiguously as the discovery of a new light spin 1 particle, in contrast with 
other experiments searching for light weakly interacting particles in rare K, ir or fj, decays |], |6|. 

From the analysis of the data from earlier experiments, constraints on the branching ratio for 
the decay of P — > 7 + X range from 10 -7 to 10 -3 depending on whether X interacts with both 
quarks and leptons or only with quarks. Since in the first case X is a short lived particle decaying 
mainly to e + e~ or vv pairs, we will only consider the second case, where X is a relatively long-lived 
particle ||. 

Direct searches for a signal from 7r° — > 7 + X decay have been performed in a few experiments 
with two different methods: i) searching for a peak in inclusive photon spectra from two-body 
7T° — > 7 + nothing decays, where "nothing" means that X is not detected because it either has a 
long life time or decays into vu pairs ( ref. H |f|); and ii) searching for a peak in the invariant 
mass spectrum of e + e~ pairs from tt° decays, which corresponds to the decay X — > e + e~ [10]. 

The best experimental limit on the decay tt° — > 7 + X was obtained recently by the Crystal 
Barrel Collaboration at CERN ||. Using pp annihilations as a source of neutral pions, they 
searched for a single peak in the inclusive photon energy spectrum. The branching ratio limit of 
(6 to 3)xl0~ 5 (90%C.L.) was obtained for 65 < M x < 125 MeV. A less stringent upper limit 
BR(tt° -> 7 + X) < (3 to 0.6) xlCT 4 was obtained for the mass region < M x < 65 MeV/c 2 . This 
result is valid for the case where X is a long-lived particle or it decays mainly to vv pairs. 

In this paper, we present a more sensitive upper limit on the branching ratio of the decay 
vr° — > 7 + X obtained by using a new method |ll| from the analysis of high energy neutrino data 
taken by the NOMAD experiment at CERN. 



2 The NOMAD Detector 

The NOMAD detector, designed to search for a neutrino oscillation signal in the CERN SPS wide- 
band neutrino beam, is described in detail in ref.[|l2|. A sketch of the NOMAD detector is shown 
in Figure |l|. It consists of a number of subdetectors most of which are located inside a 0.4 T dipole 
magnet with a volume of 7.5x3.5x3.5 m 3 . The relevant features for the present study will be briefly 
mentioned. 

The complete active target consists of 44 drift chambers (DC) mounted in 11 modules [|i"3"| . The 
target is followed by a transition radiation detector (TRD) to enhance e/n separation and by a 
lead-glass electromagnetic calorimeter (ECAL) with an upstream preshower detector (PRS). Five 
additional drift chambers are installed in the TRD region. 

Each of the 9 TRD modules [[IJJ] consists of a radiator followed by a detection plane of vertical 
straw tubes. 

The ECAL consists of 875 lead-glass Cerenkov counters of TF1-000 type arranged in a matrix 
of 35 rows by 25 columns. Each counter is about 19 radiation lengths deep. The energy resolution 
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is a/E = 0.01 + 0.032/\/i?, where E is the shower energy in GeV. A more detailed description of 
the ECAL is given elsewhere |15[| . 

The PRS is composed of two planes of proportional tubes (286 horizontal and 288 vertical tubes) 
preceded by a 9 mm thick (1.6 Xo) lead plate, providing a spatial resolution a x ^ y ~ lcm/ 'y/E, where 
E is the shower energy in GeV. The fiducial mass of this detector is about 700 kg. 

The trigger for neutrino interactions in the target is provided by two planes of scintillation 
counters Ti and T2. A veto in front of the magnet rejects upstream neutrino interactions and 
muons incident on the detector. Neutrino interactions in the PRS or ECAL are collected by a 
Ti x T2 x ECAL trigger specially designed for this purpose. The ECAL signal is obtained as the 
OR of all counter signals exceeding a threshold of about 0.8 GeV. The timing of this signal depends 
on the deposited energy. For energies above 3 GeV a time resolution of a few ns and a trigger 
efficiency of 100% are obtained. The average rate of the T\ x T2 x ECAL trigger is about 3/ 10 13 
protons on the neutrino target (p.o.t.). 

A hadronic calorimeter(HCAL) and a set of 10 drift chambers placed behind the magnet pro- 
vide an estimate of the energy of the hadronic component in the event and muon identification, 
respectively. The HCAL is an iron-scintillator sampling calorimeter consisting of 11 iron plates, 4.9 
cm thick, separated by 1.8 cm gaps in which scintillator paddles 3.6 m long 1 cm thick, and 18.3 
cm wide are installed. The HCAL active area is 3.6 m wide by 3.5 m high, and approximately 3.1 
hadronic interaction lengths deep. 



3 Method of Search 

If the decay 7r° — > 7 + X exists, one expects a flux of high energy X bosons from the SPS neutrino 
target, since 7r°'s are abundantly produced in the forward direction by 450 GeV protons either in 
the beryllium target or in the beam dump following the decay tunnel. If X is a long-lived particle, 
this flux would penetrate the downstream shielding without significant attenuation and would be 
observed in the NOMAD detector via the Primakoff effect, namely in the conversion process X —* ir° 



in the external Coulomb field of a nucleus [11] (see Figure 2). 

Because the cross section for X — > tt° conversion is proportional to Z 2 , we searched for these 
events in the lead of the preshower detector. The experimental signature of X — ► 7r° conversion is a 
single high energy ir° decaying into two photons which results in a single isolated electromagnetic 
shower in the ECAL. The corresponding ECAL cluster should be matched to the PRS cluster 
from the converted photons and must not be accompanied by a significant activity in any of the 
other NOMAD subdetectors. The occurrence of X — > tt° conversion would appear as an excess of 
neutrino-like interactions in the PRS with pure electromagnetic final states above those expected 
from Monte Carlo predictions of standard neutrino interactions. 

The expected energy spectrum of X-bosons at the NOMAD detector is shown in Figure 3 
for a mass Mx=10 MeV. The energy spectra of 7r°'s produced in the neutrino target and in the 
beam dump have been obtained with the same detailed GEANT |l(| simulation used to predict the 
neutrino flux distributions at the NOMAD detector. The X-spectrum is considerably harder than 
that of neutrinos which is also shown for comparison. 



4 Data Sample and Event Selection 

This analysis is based on the data taken during the first half of the neutrino run in 1995, in which 
the NOMAD target was only partially installed and consisted of four drift chamber modules placed 
upstream of the TRD detector. The integrated number of protons delivered to the neutrino target 
during this period was about 1.45 xlO 18 . 



4 



Monte Carlo simulations use the LEPTO 6.1 event generator [JD]] supplemented by event gen- 
erators for resonance, quasi-elastic and coherent neutrino processes and a full detector simulation 



based on GEANT[16j. These simulations together with direct measurements of subdetector occu- 
pancy during the neutrino spill with a random trigger [^] show that v and X — ► 7r° events in the 
PRS are accompanied by no significant activity in the DC and TRD. For this reason a simple cut on 
the number of hits iV^f, N^ (TRD) and N^ D in the DC, in the DC located in the TRD region, 
and in the TRD, respectively, were used as a veto to suppress events from neutrino interactions 
which occurred in these regions or in the magnet coils. However, neutrino events, mainly with a 
neutral final state, which occurred in the near upstream PRS region can also pass these veto cuts 
and can be taken as neutrino events in the PRS. For this reason the upstream region was also 
included in the simulations and surviving events were considered as neutrino events in the PRS. 

The selection criteria for X — ► ir° events are based on a full Monte Carlo simulation of X — > tt° 
conversions in the NOMAD detector and rely on their properties as mentioned in Section 3, namely: 

• the quality of the match between the ECAL shower and the corresponding PRS cluster; 

• the total energy and the shower shape in the ECAL; 

• the amount of energy deposition in the HCAL. 

Candidate events were identified by the following simultaneous requirements: 



• DC x TRD (no activity in the DC or TRD ): JVjgf < 10, iVjJJ ; ' < 4, N^f D < 2. 

• PRS x ECAL(> 8 GeV): isolated PRS cluster matched to isolated ECAL cluster in both X 
and Y planes. 

At most two ECAL clusters in the ECAL were allowed. The energy of the most energetic 
cluster had to be Eecal > 8 GeV and the second cluster had to be less than 0.3 GeV in order 
to reject pile-up events and ECAL noise. The shape of the most energetic cluster was fitted 
to the shape expected from an electromagnetic shower and the \ 2 °f the fit was required 
to be less than 20 J18[ |. The differences AX and AY between the X and Y coordinates 
of the cluster in the PRS and the corresponding cluster in the ECAL were required to be 
\AX\ < 2 cm, \AY\ < 1.5 cm. 

These conditions were used to identify isolated electromagnetic showers in the ECAL that 
originated from the conversion of photons from 7r°'s produced in the preshower. 



• HCAL (no HCAL activity): the total visible energy was required to be less than the HCAL 
noise threshold Ehcal < 0.4 GeV. 

This cut serves to identify electromagnetic energy in the ECAL. The reliability of the HCAL 
veto to suppress conventional neutrino events with an energy leakage to the HCAL detector 
was checked with straight through muons and with a selected sample of Un charged current 
(CC) events in the PRS. 



• MUON : no track(s) in the muon chambers matched to the PRS cluster. 

After applying these cuts to the initial sample of 4.83 x 10 5 events recorded with the ECAL 
trigger we found 20 candidate events for X — > tt° conversion. The amount of background in this 
sample from standard and v e interactions was evaluated using the Monte Carlo (see Section 5). 
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Applying the DC x TRD x PRS x ECAL(> IGeV) x MUON selection criteria to the same 
initial sample we find 3691 events from v^CC interactions in the PRS. Here, MUON denotes 
single muon tracks extrapolated back from the muon chambers to the PRS and matched to the 
PRS cluster. The spectra of energy deposited in the ECAL for these selected v^CC events and for 
simulated events which pass the same reconstruction program and selection cuts are found to be 
in agreement. The T\ x T2 x ECAL trigger efficiency was obtained from a Monte Carlo simulation 
and was found to be 71% and 97% for v^CC and X — > tt° events, respectively. Taking the overall 
selection efficiency into account, the total number of v^CC events interacting in the PRS was found 
to be iV„ CK7=2.29 x 10 4 . The number of v^CC events was converted to the number of protons on 
target N po t using the simulation of the neutrino flux at the NOMAD detector and the known cross 
section value for v^CC interactions, giving Np t = 1-36 x 10 18 . This number agrees within 10% with 
the value of N pot measured by the beam monitors and was used further for normalisation. The main 
uncertainty results from the contribution of events from neutrino interactions in the downstream 
TRD region which pass the DC and TRD cuts and from backscattering in events occurring in the 
PRS. By varying the DC and TRD cuts, it was found that the systematic error in the number of 
ffj.CC events in the PRS resulting from backscattering is of the order of 10%. 

5 Background Events 

The main background to X — > 7r° conversions is expected from the following neutrino processes 
occurring in the PRS or upstream PRS region with a significant electromagnetic component in the 
final state and with no significant energy deposition in the HCAL : 

• v^CC interactions classified as muonless because the muon was not detected; 

• inclusive 7r° production from neutral current (NC) interactions; 

• coherent and diffractive 7r° production; 

• quasi-elastic v e scattering. 

• v e CC interactions; 

Backgrounds from the trident events, v e + Pb — ► v e + e + + e~ + Pb and from pile-up events 
were found to be negligible. 

To evaluate the amount of background in the data sample, simulated events were processed 
through the same reconstruction program and selection criteria that were used for real neutrino 
data. Then, all background distributions were summed up, taking into account the corresponding 
normalisation factors. These factors were calculated from beam composition and cross sections of 
the well-known processes listed above. The total number of events interacting in the PRS and the 
number of expected candidate events after applying the selection criteria are given in Table 1 for 
each background process. 

The total background in the data sample was estimated to be 18.1±2.8 events, where statistical 
and systematic errors are added in quadrature. The fraction of neutrino interactions in the PRS 
which satisfy all cuts is 6 x 10~ 4 . 

6 Result 

Figure 4 shows the combined background and candidate event energy spectra in the ECAL. The 
agreement between data and Monte Carlo is good. The overall efficiency for X — > ir° conversion 
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Table 1: Overall background estimate 



item 


VJCC 




Coh.^ u 


u e - QE 


v e CC 


Total number of 
interactions 
in the PRS 


22886 


6866 


59 


6.2 


320 


Number of 
expected 
candidate events 


1.7±1.0 


6.6±1.3 


1.5±0.4 


2.1±0.2 


6.2±1.2 



detection was found to be ~24%. The inefficiency is mostly due to the requirement that at least 
one photon from tt° decay convert in the PRS. 

By subtracting the number of expected background events from the number of candidate events 
we obtain Nx^ n o = 1.9±5.3 events showing no excess of X — ► 7r° conversion-like events, and, hence, 
no indication for the existence of this process. The 90% C.L. upper limit for the branching ratio 
Br(TT° — ► 7 + X) was calculated by using the following relation: 



Nx"^ > [Br(ir v -> 7 + X)] z x J ^—a (M x , E x )dE x x e sel x pi— (1) 

where JVjj?t° (=8.7 events) is the 90% C.L. upper limit for the expected number of sig- 
nal events, e se / is the selection efficiency, which was found to be practically independent of Ex, 
fo(Mx,Ex, Np t) and ao(Mx, Ex) are the flux of X bosons for the given number of protons Npot 
on the neutrino target and the cross section for X — » ir° conversion on lead, respectively, calculated 
for Br(ir° — > 7 + X)=l. In Eq.(l), Mx and Ex are the mass and energy of the X boson, respec- 
tively. The cross section oq(Mx,Ex) is given in ref. [11|. The total X— flux per proton on target 



calculated as a function of the X boson mass is shown in Figure 5. We note that J Br(vr° -> 7 + X) 
appears twice in the formula for Nx^^o, through the X boson flux from the target and through 
the Primakoff mechanism. 

The 90 % C.L. branching ratio limit curve is shown in Figure 6 together with the result of 
ref. . For the mass region 0< Mx < 120 MeV/c 2 the limit is 

Br(TT° -> 7 + X) < (3.3 to 1.9) x 10~ 5 (2) 

Varying the cut on the ECAL energy deposition in the range from 5 to 50 GeV does not change the 
limit substantially, while above 50 GeV it becomes worse. The limit is valid for an X boson lifetime 
tx > W~ 9 Mx[MeV/c 2 ] s. For the mass region 0< Mx < 60 MeV/c 2 , the limit is approximately 
a factor 10 to 5 better than the best previously published limit obtained by the Crystal Barrel 
collaboration Q. 

Our result can also be used to constrain the magnitude of the coupling of a hypothetical gauge 
X boson to quarks || [ll]] 



g 2 < 2 x 10" 7 (1 - ^f)" 3 (3) 
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The attenuation of the X- flux due to X interactions with matter was found to be negligible, 
since for Br(n° — > X7) < 10~ 4 the X boson mean free path in iron is > 300 km, as compared 
with the iron and earth shielding total length of 0.4 km used in our beam. Furthermore the decay 
length for Ex > 10 GeV was estimated to be longer than ~10 km [O]. We note that limits on the 
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branching ratio Br(rj, rj — > 7 + X) could also be obtained from our data if the cross-sections for 
77, rf productions in the forward direction in p — Be collisions at 450 GeV were known. 
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Figure 1: Side view of the NOMAD detector. 
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Figure 3: Combined energy spectrum of X bosons with mass Mx=10 Me V from the SPS 
neutrino target and from the beam dump region at the NOMAD detector. The neutrino 
energy spectrum (hatched area) is also shown for comparison. The two spectra are arbitrarily 
normalized. 
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Figure 4: ECAL energy distributions for candidate events (shaded area) and for combined 
background events. 
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Figure 5: Total X boson flux at the NOMAD detector versus mass Mx ( Ex > 8 GeV) 
calculated for Br(n° — > 7 + X)=l. The curve is a polynomial fit to the points. 



14 



-3 




20 40 60 80 100 120 

M x (MeV/c 2 ) 



Figure 6: 90% confidence level upper limit on branching Br(n° — » 7 + X) versus Mx- The 
limit is valid for a lifetime of the X boson tx > 1CT 9 Mx[MeV/c 2 ] s. 
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